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Abstract: We derive in a straightforward way the null controllability of a 1-D heat equation 
with boundary control. We use the so-called flatness approach, which consists in parameterizing 
the solution and the control by the derivatives of a "flat output". This provides an explicit control 
law achieving the exact steering to zero. We also give accurate error estimates when the various 
series involved are replaced by their partial sums, which is paramount for an actual numerical 
scheme. Numerical experiments demonstrate the relevance of the approach. 
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1. INTRODUCTION 

The controllability of the heat equation was first consid- 
ered in the 1-D case, Fattorini and Russell (1971); Jones Jr. 
(1977); Liftman (1978)), and very precise results were 
obtained by the classical moment approach. Next using 
Carleman estimates and duality arguments the null con- 
trollability was proved in Fursikov and Imanuvilov (1996); 
Lebeau and Robbiano (1995) for any bounded domain 
in WL N , any control time T, and any control region. This 
Carleman approach proves very efficient also with semilin- 
ear parabolic equations, Fursikov and Imanuvilov (1996). 

By contrast the numerical control of the heat equation (or 
of parabolic equations) is in its early stage, see e.g. Munch 
and Zuazua (2010); Boyer et al. (2011); Micu and Zuazua 
(2011). A natural candidate for the control input is the 



control of 



ral L — 



norm, w 



hich may be obtained 



as a trace of the solution of the (backward) adjoint 
problem whose terminal state is the minimizer of a suitable 
quadratic cost. Unfortunately its computation is a hard 
task Micu and Zuazua (2011); indeed the terminal state 
of the adjoint problem associated with some regular initial 
state of the control problem may be highly irregular, which 
leads to severe troubles in the numerical computation of 
the control function. 

All the above results rely on some observability inequalities 
for the adjoint system. A direct approach which does not 
involve the adjoint problem was proposed in Jones Jr. 
(1977); Littman (1978); Lin Guo and Liftman (1995). 
In Jones Jr. (1977) a fundamental solution for the heat 
equation with compact support in time was introduced and 
used to prove null controllability. The results in Jones Jr. 
(1977); Rosier (2002) can be used to derive control results 
on a bounded interval with one boundary control in 
some Gevrey class. An extension of those results to the 
semilinear heat equation in ID was obtained in Lin Guo 



and Littman (1995) in a more explicit way through the 
resolution of an ill-posed problem with data of Gevrey 
order 2 in t. 

In this paper we derive in a straightforward way the null 
controllability of the 1-D heat equation 

0t(t,x)-9 xx (t,x) = 0, (t,x) e (0,T) x (0,1) (1) 

x (t,O) = O, te(0,T) (2) 

x (t,l) = «(*)> te(0,T) (3) 
with initial condition 

0(O,x) = 9 o (x), xe(0,l). 

This system describes the dynamics of the temperature 
in an insulated metal rod where the control u is the heat 
flux at one end. More precisely given any final time T > 
and any initial state 0o € L 2 (0, 1) we provide an explicit 
control input u G L 2 (0,T) such that the state reached at 
time T is zero, i.e. 

0(T,x)=O, see (0,1). 

We use the so-called flatness approach, Fliess et al. (1995), 
which consists in parameterizing the solution and the 
control u by the derivatives of a "flat output" y (sec- 
tion 2); this notion was initially introduced for finite- 
dimensional (nonlinear) systems, and later extended to in 
particular parabolic PDEs, Laroche et al. (2000); Lynch 
and Rudolph (2002); Meurer and Zeitz (2008); Meurer 
(2011). Choosing a suitable trajectory for this flat output y 
then yields an explicit series for a control achieving the 
exact steering to zero (section 3) . This generalizes Laroche 
et al. (2000), where only approximate controllability was 
achieved through a similar construction. We then give 
accurate error estimates when the various series involved 
are replaced by their partial sums, which is paramount 
for an actual numerical scheme (section 4). Numerical 
experiments demonstrate the relevance of the approach 
(section 5). 



In the sequel we will consider series with infinitely many 
derivatives of some functions. The notion of Gevrey order 
is a way of estimating the growth of these derivatives: we 
say that a function y G C°° ([0, T]) is Gevrey of order s > 
on [0, T] if there exist positive constants M, R such that 



y (p) {t) 



< M 



p\° 



V£ G [0,T], Vp > 0. 



More generally if K C R N is a compact set and y is a 
function of class C°° on X (i.e. y is the restriction to K 
of a function of class C°° on some open neighbourhood £2 
of -ftT), we say y is Gevrey of order s\ in x\, s 2 in x 2 ,. ■ ■ ,sm 
in xn on K if there exist positive constants M,R\,..., Rn 
such that 

Yl N (n-l) Si 

■ ■ ■ d™y(x)\ < ^ u p- ■ g V P g N w . 
1 li=i H i 

By definition, a Gevrey function of order s is also of order 
r for r > s. Gevrey functions of order 1 are analytic 
(entire if s < 1). Gevrey functions of order s > 1 have 
a divergent Taylor expansion; the larger s, the "more 
divergent" the Taylor expansion. Important properties of 
analytic functions generalize to Gevrey functions of order 
s > 1: the scaling, addition, multiplication and derivation 
of Gevrey functions of order s > 1 is of order s, see Ramis 
(1978); Rudin (1987). But contrary to analytic functions, 
functions of order s > 1 may be constant on an open set 
without being constant everywhere. For example the "step 
function" 

( 1 if £ < 

if £ > 1 

e -(i-t)- fe 

, e -(i-«)-» +e -«-» iftG]0 < 1[ < 

1 is Gevrey of order s on [0, 1] (and 
(0) = 1, 0,(1) = and #(0) = 



where k = (s — 1) 
in fact on R); notice 



(1) = for all i > 1. 



In conjunction with growth estimates we will repeatedly 
use Stirling's formula n\ <~ (nje) n \phm. 

2. THE HEAT EQUATION IS FLAT 

We claim the system (l)-(3) is "flat" with y(t) := 6(0, t) 
as a flat output, which means there is (in appropriate 
spaces of smooth functions) a 1 — 1 correspondence between 
arbitrary functions 1 1-> y(t) and solutions of (l)-(3). 

We first seek a formal solution in the form 

x l 



6{t,x) == -Tj-OiW 



i>0 



where the a^s are functions yet to define. Plugging this 
expression into (1) yields 



^7r[ a *+ 2_a ^ =0 ' 



hence a i+2 — a- for all i > 0. On the other hand y(t) = 
9(0, t) = ao(t), and (2) implies a\(t) — 0. As a consequence 
a 2i = and d2i+i = for all i > 0. The formal solution 
thus reads 

(2z)! ; 



(4) 



while the formal control is given by 



i>i v ; 



(5) 



We now give a meaning to this formal solution by restrict- 
ing 1 1-> y(t) to be Gevrey of order s G [0, 2). 
Proposition 1. Let s G [0,2), — oo < t\ < t 2 < oo, and 
y G C°°([ti, i 2 ]) satisfying for some constants M, R > 



y [l) (t) 



< M- 



R l 



Mi > 0, Vt G [ti,t 2 \. 



(6) 



Then the function 9 defined by (4) is Gevrey of order s in t 
and s/2 in x on [ti,t 2 ] x [0, 1]; hence the control u defined 
by (5) is also Gevrey of order s on [ii,i2]- 



Proof. We must prove the formal series 

Jli-n 



d?d n j{t, X ) = Y J 



(2i-n)\ 



(t+n) 



(*) 



(7) 



is uniformly convergent on [£i , £2] x [0, 1] with growth 
estimates of the form 

ITO(M)|<C^|- (8) 

By (6), we have for all (t,x) G [£1,^2] x [0, 1] 



(2i-n)\ 



(t+m) 



(*) 



< 



< 



M (i + m)! s 

i? i+m (2« - n)! 

M (2 l+m i\m\) s 
R i+m (2i - nj\ 



< 



M 2 si (2- 2i V^i(2i)\y 



R i+r 

<M- 



(2i - n)\ 



Q — sm 



(iri)' 



— n' 2 - 
i? l 1 (2i-n)! 1 -§ ' RT' 

where we have set R\ — 2~ S R; we have used Stirling's 
formula for (2i)\ and twice (i + j)\ < 2' 1+ H\j\. Since 

< 00 the series in (7) are uniformly 



E 

2i>n 



i>0 



A^2i>n K i( 2i _„)! 1 -i 

convergent for all m, n > 0, hence 9 G C°°([£i, t 2 ) x [0, 1]). 
Finally, since 

where R 2 G (0, \/R\) and C > is some constant 
independent of n, we have the desired estimates (8). □ 

3. NULL CONTROLLABILITY 

In this section we derive an explicit control steering the 
system from any initial state # G L 2 (0,1) at time to 
the final state at time T > 0. Two ideas are involved: 
on the one hand thanks to the flatness property it is easy 
to find a control achieving the steering to zero starting 
from a certain set of initial conditions (lemma 2); on the 
other hand thanks to the regularizing property of the heat 
equation this set is reached from any 9 G L 2 (0, 1) when 
applying first a zero control for some time (lemma 3). 

Lemma 2. Let (j/j)j>o be a sequence of real numbers such 
that for some constants M,R> 

\Vi\<M± V z>0. (9) 



Then the function defined on [ti, t\ + R'], R! < R, by 



i>0 



is Gevrey of order s > 1 on [f i, t\ + R'] and satisfies for all 
i > 

y (l Hh)=y l (10) 

y^(t 1 +R')=0. (11) 



Moreover the control defined on \t\,t\ + R'] by 



i>l 



(12) 



is also Gevrey of order s on + R'] and steers the 

system from the initial state ^2 i>0 ^* (fi)T a * *i *° t ne 
final state at time t\ + R' . 

Proof. Let y(z) := X)»>of*fr- The g rowtn property (9) 
implies y is analytic on the disc {z G C; \z\ < R}, the 
convergence being moreover uniform for \z\ < R! < R. 
Therefore y is Gevrey of order 1, hence of order s > 1, 

on [<i,ti+_R']. On the other hand (j>{t) := 4>s(jj^j is also 

Gevrey of order s on \t\ , ti + R'] , hence so is the product y 
of y and (j). The boundary values (10)-(11) follow at once 
from the definition of <fr s . 

The control u in (12) achieves the steering to zero; indeed 



i>0 



as well as u, is by proposition 1 Gevrey of order s in t and 
s/2 in x, and obviously satisfies 6(t\ + R' , x) = 0. □ 

Lemma 3. Let O € L 2 (0, 1) and t > 0. Consider the final 
state 9 T (x) := 9(t,x) reached when applying the control 
u(t) := 0, t € [0, r], starting from the initial state 9 - 

Then # T is analytic in C and can be expanded as 

~1i 



i>0 



with 



where C is some positive constant depending only on 9 . 
Proof. Decompose 9 as the Fourier series of cosines 
6o{ x ) = E c n v / 2cos(n7rx) 



n>0 



where the convergence holds in L 2 (0, 1) and 



2M 



n>l J ° 



x)\ 2 dx < oo. 



The solution starting from 9 then reads 

9(t,x) = ^c„e-™ 27r2t v / 2cos(n 7 rx) 



(13) 



n>0 



and in particular 

^0*0 = ^c n e- n2 ^ T V2cos{m:x). 



The series for T is analytic in C since for all |x| < r 

c n e- n27l2T V2cos(mrx) < C\ (sup\c k \) e -™ 2 ^+™?"- 

\fe>o / 

where C\ is some positive constant; this ensures the 
uniform convergence of the series in every open disk of 
radius r > 0. 



Moreover 

T (x) = y2^ C „e- n2 - 2 ^(-l) f 



n>0 
„2i 



i>0 



(rmx) 2 ' 1 



i>0 



C n 6 



>7T) 2 



n>0 



The change in the order of summation will be justified once 
we have proved that yt, i > 0, is absolutely convergent and 

„2i 



El* 



i>0 



(2i)! 



< oo, Vx > 0. 



For i > let h t (x) := e^" 1 ^ 2 (the) 2 * and iV 4 := (^) 

i 

The map hi is increasing on [0, ( 7r i J :) 2 ] and decreasing on 
[(^if) 5 ,+oo) hence 



n>0 



V/ii(n)</ hi(x)dx + hi(Ni) 



+ /i i (iV i + l) + 



hi(x)dx 



N,. 



< c 2 



T*%A ./(| 



hi(x)dx; 



Ci is some positive constant and we have used Stirling's 
formula. On the other hand integrating by parts yields 



/ hi(x)dx = / hi_i(x)dx 

Jo 2t J 

(2i — 1) - - ■ 3 - 1 



(2r)< 
(2t)! 1 
24!(2t)* 7r^7 



dx 



e x dx 



<c 3 — =, 

T l \l IT 

where C 3 is some positive constant and we have again used 
Stirling's formula. As a consequence 

|y i |<\/2sup|c„|^/i i (n)<c(l + ^)-^ (14) 

n>0 



n>0 



(2t 



y 



where C is some positive constant. Finally 
v^, , x 21 / 1 n — i! 2 

2^ If' 

j>0 



00 



since 



Vi 4i t 



□ 



n>0 



With the two previous lemma at hand we can now state 
our main controllability result. 



Theorem 4. Let E L 2 (0, 1) and T > 0. Pick any 
t e (0, T) and s E (1,2). Then there exists a function y 
Gevrey of order s on [r, T] such that the control 



u(t) := 





E 



i/ (<) (*) 



t (2i-l)! 



if < i < t 
if t < i < T. 



steers the system from the initial state 9q a t time to the 
final state at time T. 

Moreover u is Gevrey of order s on [0,T]; t i-> •) is in 
C([0, T],L 2 (0, 1)); is Gevrey of order s in t and s/2 in ir 
on [e,T] x [0,1] for all e € (0,T). 

Proof. By lemma 3 the state reached at time r reads 



i>0 J/» (2»)l 



with the sequence (y. 



growth property of lemma 2 with M := 1- 
Hence the desired function is given by 



satisfying the 
-f= and R := r. 



(t-T)' 



i>0 



on [t, t + i?'], where R' < t and i?' < T - r; and 
by y(t) := on [r + i?',T]. Moreover y and u are 
Gevrey of order s on [r, T], and by construction 9(t,x) = 

0{t+,x) = E^o^tSr and u(t) - u(t+) = 0. By 
Proposition 1 the solution on [r, T] x [0,1] is well- 
defined and Gevrey of order s in t and s/2 in x hence 
1 1 y 9(t,-) E CftO.Tj'jC^O,!])) and u E C([0,T]). On the 
other hand it is easily seen that 9 is Gevrey of order 1 in t 
and 1/2 in x on [e,r] x [0, 1] for all e E (0,r). 

Thus the solution # is Gevrey of order 1 in t and 1 /2 in x 
on [e,r] x [0, 1] while it is Gevrey of order s in t and s/2 
in x on [r, T] x [0, 1]. To prove 9 is Gevrey of order s in t 
and s/2 in .t on [e,T] x [0, 1] it is then sufficient to check 
d?9(T, x) = d^9{T+,x) for k > and a; G [0, 1]. But 



i>0 



E 

i>0 



(2i)!' 



(»+fc) 



«(*) = E (2731)1 
»>1 v y 



i>0 



(2i)l 



2/fe- 



fc>0 



{t-rf 
k\ 



(15) 
(16) 
(17) 



y k = yfife c n e- n2 ^n 2k ) (-7r 2 ) fc ; (18) 



n>0 



moreover y hence u and are identically zero on [t+R' , T] . 
The aim of this section is to show that the partial sums 



Ki<i 



(2i - 1) 



a; 



2i 



9(4,*):= £ y«(i)- 

0<i<i 



(t-r)» 



k\ 



(19) 
(20) 
(21) 



0<fe<fc 

y- k :=V2( £ ^-" 2 ^n 2fe )(-7r 2 ) fc . (22) 

0<n<n 

for given i,k,n € N provide very good approximations of 
the above series, and to give explicit error estimates. 

Theorem 5. There exist positive constants C,Ci,C2,Cs 
such that for all 9 E L 2 , i, k,n E N, and t E [t, T] 

\\9(t)-6(t)\\ LOC < C (e~ c ^ nni + e~ c ^ + e~ c ^ \\9 \\ L , 



Proof. First notice that for (t, x) E [r, T] x [0, 1] 

\9(t,x)-9(t,x) \ < A1 + A2 + A3, 

where 

A,:= E^H')E^ (/ ' " 



fc] 



0<i<i 



k>k 



(2»)! 



A, := 



E E 



0<i<i 



(2i)! 



0<fe<fc 



i>0 v ' \n 



V2(J2c n e- n ^n 2k )(-n 2 y 



: (t-r) h 



>0 



n>n 



k\ 



(2i)\ 



= J2 c n(-n 2 Tr 2 ) k e- n27r2T V2cos(mTx) 

n>0 

= d k 9(r,x). 

As a consequence u is also Gevrey of order s on [0, T]. □ 
4. NUMERICAL ESTIMATES 



Summarizing the previous section the control u and solu- 
tion 9 on [t, t + R'] are given by the infinite series 



By lemma 2 y is Gevrey of order s on [r, T] with 
some Mi ||#o|Il2 , i?i > hence 



Ar<E 



y (<) W| 

(2i)! 



<M(||^o|| L 2E 

<^iii^oii i2 E 



(2i)\R\ 



y/i /Z\( s -2)i 



(;) 



(4i?x)^ 2 " ve 

i £ T i e ( 2 - s )( 1 - lnl ) 4 

(4Ri)' : 



where we have used Stirling's formula. Pick C\ < 2 — s 
and a E [C\ , 2 — s) . Then for i > i 



(4iii)- i iV e (2-«)(l-ln^ < Kie -ai(lni-l) ^ 

where K\ = Ki(s,a,R\). But 

/•OO 

^ e -ai(lni-l) < / e - CT x(lnx-l) da , 

/•OO 

< / e- CTX da; 

•/i(lni-l) 

so that eventually 

Ai <K?\\0 o \\v, e~ c ^\ 

For A2 we first notice that for t e {z e C; \z — t\ < p'r}, 
where p' satisfies R' /t < p' < 1, 

IEf*T^I<Ei»* lP ' 



k>k 



k\ i-^_<^< k\ 

k>k 



k>k 
/fc+1 



<<||fl || i2 P' 

By the Cauchy estimates we thus have for r < t < t + i?' 



k>k 



(t ~ T)fc l < ILf'llfl II n' fe+1 * ! 

< Ml ||«0|| L 2 P — : 



jfe! 



hence 



<M 2 ll0oll i2 p' 



/fe+1 « 



fe>fe 

It follows that 



A 2 < J2 M 2\\ o\\ L 2P 



fe+1 «! 



0<i<i 



(2i)\Rl 



< K 2 \\e \\ L 2 e 



-C 2 k 



where < C 2 < In . 

To estimate A 3 first notice that for a > 

E e 



""" < / e~ ax dx = 

Jn Jn 2 



n > n 



2^ 



2 cm 



Pick p' and p" with Rt < p' < p" < 1. Then 

e _ ' c fc fe 
(tt 2 p"t)* 



E c ' 



^w^^-^E c - B * (lv ' )T 

fcl e -7r 2 (l-p")rn 2 

yK A p"T) K n 

where we have used first that x — > e~ x x k is maximum for 
x = k, and then that e~ k k k < Ck\ by Stirling's formula. 
Hence for z e {z e C; |z — r| < p'r} 



0<fe<fe ™>™ 

<^ 3 KII L 2 E 



fc (z-T) fc 



-7r 2 (l-p")Tn 2 . / . fc 



<^3 INL 2 



n 

0<fe<fe 

-7r 2 ( 1 — p")rn 2 



(50 



Setting C* 3 := tt 2 (1 - p")r < tt 2 (t - R'), this yields for 
t < t < t + R' 



9i[ E V2(^ C „ e - 2 ^)(-. 2 ) 

0<fe<fe ™>" 

2 i! 



fe (*-T) fc - 
fc! ■ 



<M 3 ||0 o || i2 e 



-c 3 ™ 2 



and 



Q<k<k ™>™ 

<M^||0 o || i2e - c - 2 iP 



fc! 



We then conclude 
A 3 < J2 KPoWl^' 



R 9 



^ 2 ^1_.<M^||0 O || L2 



0<i<i 



(2i)!i£ 



Collecting the inequalities for Ai, A 2 , A3 eventually gives 
the statement of the theorem. □ 

Let 9 denote the solution of (l)-(3) where the "exact" 
control u (15) is replaced by the truncated control u (19), 
still starting from the initial condition 9q. Notice 9 is the 
"exact" solution obtained when applying the truncated 
control, while 9 is the truncated solution obtained when 
applying the "exact" control. Since in "real life" only the 
truncated control can be actually applied it is important 
to know how well 9 approximates 9, the "exact" solution 
when applying the "exact" control. It turns that 9 satisfies 
the same relation as 9 in Theorem 5; the proof is omitted 
for lack of space but follows from the proof of Theorem 5. 

Corollary 6. With the same notations as in Theorem 5, 



9(t)-9(t) 



< C (e- Cl + e-° 2 k + e- Cs n ' 



) Pol 



L 2 



5. NUMERICAL EXPERIMENTS 



We hare conducted a number of numerical experiments 
to demonstrate the relevance of our approach, and to 
investigate the influence of the main parameters. We have 
focused on the control effort, which is probably the most 
important quantity. The results are summarized in the two 
tables below, where ||u|| i2 (top table) and ||u|| iDO (bottom 
table) are given for different values of the parameters s 
and R' in (21). Figures 1 and 2 give moreover the complete 
temperature and control evolution for the case R! = 0.2 
and s = 1.6. For all experiments the "regularization time" 
t is 0.3, and the initial condition 9q is a step function with 
9 (x) = -1 on [0,1/2) and 9 (x) = 1 on (1/2,1) hence 

Fourier coefficients c 2p = and c 2p+ i = ^p+i ^ or 
p > 0; i,k,n are "large enough" for a good accuracy of the 
truncated series(19) -(22). 



s,R' 


0.15 


0.20 


0.25 


0.30 


1.5 


693 


63.3 


12.7 


3.82 


1.6 


35.3 


6.41 


1.95 


0.78 


1.7 


7.49 


1.95 


0.74 


0.34 


1.8 


5.53 


1.24 


0.48 


0.23 


1.9 


5.71 


1.29 


0.47 


0.22 



s,R' 


0.15 


0.20 


0.25 


0.30 


1.5 


3666 


330 


55.2 


18.1 


1.6 


118 


23.6 


7.17 


2.76 


1.7 


18.6 


4.78 


1.73 


0.76 


1.8 


36.4 


4.59 


1.44 


0.65 


1.9 


47.8 


9.66 


2.13 


0.82 



A few conclusions can be drawn from these experiments: 

• when R' gets small there seems to be an "optimal" 
value of s which minimizes the control effort (the 
values are different for ||u|| i2 and 

• when s it too small or too large the derivatives y( % ' 
tend to "crowd" near the extremities of [t, t+R'] with 
large amplitudes 

• the smaller R' the more noticeable this phenomenon, 
with of course an increased control effort. 

An interesting question is the tradeoff between r and R' to 
reach the zero state at time T := t + R' with the smallest 
control effort: longer regularization r or longer duration R' 
of the active control? With the present construction of the 
function y (18) we have by design R' < r, which is proba- 
bly an important restriction. Other, but more complicated, 
constructions without this limitation are possible and will 
be studied in the future. 




time t 



Fig. 1. 9(t,x) for R' = 0.2 and s = 1.6. 
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